Concurrence of Arbitrary Dimensional Bipartite Quantum States 
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We derive an analytical lower bound for the concurrence of a bipartite quantum state in arbitrary 
dimension. A functional relation is established relating concurrence, the Peres-Horodecki criterion 
and the realignment criterion. We demonstrate that our bound is exact for some mixed quantum 
states. The significance of our method is illustrated by giving a quantitative evaluation of entangle- 
ment for many bound entangled states, some of which fail to be identified by the usual concurrence 
estimation method. 
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Entanglement is a striking feature of quantum sys- 
tems and is the key physical resource to realize quan- 
tum information tasks such as quantum cryptography, 
quantum teleportation and quantum computation 
which cannot be accounted for by classical physics. This 
has provided a strong motivation for the study of de- 
tection and quantification of entanglement in an oper- 
ational way. Despite of a great deal of effort in past 
years [1 11 11 11 11 0, S i, [ll IH III El for the mo- 
ment only partial solutions are known for generic mixed 
states. As for quantitative measures of entanglement, 
there is an elegant formula for 2 qubits in terms of con- 
currence, which is derived analytically by Wootters in 
Ref. 0]. This quantity has recently been shown to play 
an essential role in describing quantum phase transition 
in various interacting quantum many-body systems 0| 
and may affect macroscopic properties of solids signifi- 
cantly |l5| . Furthermore, value of concurrence will pro- 
vide an estimation |l3l| for the entanglement of forma- 
tion (EOF) [3, which quantifies the required minimally 
physical resources to prepare a quantum state. It is thus 
very important to have a precise quantitative picture of 
entanglement in order to get a better insight into the 
corresponding physical systems. 

However, calculation of the concurrence is a formidable 
task as the Hilbert space dimension is increasing, like in 
the case of two parts in a real solid-state system con- 
sidered for quantum computation. Good algorithms and 
progresses have been obtained concerning lower bounds 
for qubit-qudit system ll Cl |l ll| and for bipartite systems 
in arbitrary dimension p|, ll.'^ . Considerable progress is 
made in [13 to give a purely algebraic lower bound. Nev- 
ertheless, an optimized bound generally involves numer- 
ical optimization over a large number of free parameters 
in a level (at least m{m— l)n{n— l)/4 for a m®n bipar- 
tite system, where m, n are Hilbert space dimension for 
two subsystems respectively [1,0,E.3J)- This leads to a 
computationally untractable problem for realistic system 
with a higher dimension. In addition, these methods for 
evaluating concurrence can not detect reliably arbitrary 
entangled states even if one applies all known optimiza- 
tion methods (l,*^ . 



Our aim in this work is to improve this situation dra- 
matically by giving an analytical lower bound for concur- 
rence of any mixed bipartite quantum state. We find an 
essential quantitative relation among this measure and 
available strong separability criteria. A functional rela- 
tion is explicitly derived to give a tightly lower bound for 
the concurrence. It is shown to be exact for some special 
class of states. Our method is further demonstrated to be 
better than the regular method for concurrence optimiza- 
tion, in the sense that it can detect and give an evaluation 
of entanglement for many bound entangled states (BES) 
which cannot be identified by the latter. This also com- 
plements a number of existing methods involving numer- 
ical optimization and provides a computational method 
to estimate manifestly the actual value of concurrence for 
any bipartite quantum state. 

We start with a generalized definition of concur- 
rence for a pure state |'0) in the tensor product TiA^Ti-B 
of two (finite dimensional) Hilbert spaces Hat'Hb for 2 
systems A,B. The concurrence is defined by C{\ip)) — 
^2(1 — Trp^), where the reduced density matrix pA is 
obtained by tracing over the subsystem B. The concur- 
rence is then extended to mixed states p by the convex 
roof. 
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for all possible ensemble realizations p = '^^Pi\'4'i){''Pi\: 
where pi > and X^iPi — 1- For any pure product 
state Itp), C{\ip)) vanishes according to the definition. 
Consequently, a state p is separable if and only if C (p) — 
and hence can be represented as a convex combination 
of product states as p = J2iPiPt ® pf where pf and pf 
are pure state density matrices of the subsystems A and 
B, respectively [lij . 

The key point of our idea is to relate directly the con- 
currence and the Peres-Horodecki criterion of positivity 
under partial transpose (PPT criterion) 0, 01 and the re- 
alignment criterion HQ by means of Schmidt coefficients 
of a pure state. Let us firstly consider the concurrence 
for a pure state. CdV'V) is invariant under a local unitary 
transformation (LU) 0, 0| . Without loss of generality, 
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we suppose that a pure m ® n (m < n) quantum state 
has the standard Schmidt form 

l^> =I]V7^|a.M, (2) 

i 

where ^JJTi (i = 1, . . .m) are the Schmidt coefhcients, \ai) 
and \bi) are orthonormal basis in Ha and Hb, respec- 
tively. It is evident that the reduced density matrices pA 
and pb have the same eigenvahies of /i^. It follows 

c^im^2{l-J2^^)=^J2^'^^^,, (3) 

i i<j 

which varies smoothly from 0, for pure product states, to 
2(to — l)/m for maximally entangled pure states. 

In order to derive a quantitative connection with the 
PPT criterion and the realignment criterion, we recall 
some details of the two criteria. Peres made firstly an 
important step forward for separability criterion in 2] 
by showing that p'^^ > should be satisfied for a sep- 
arable state, where p'^^ stands for a partial transpose 
with respect to the subsystem A. p^^ > is fur- 
ther shown by Horodecki et al. Q to be sufficient for 
2x2 and 2x3 bipartite syst ems. llp'^'^H is LU in- 
variant as shown in Refs.0, 0| where || • || stands for 
the trace norm defined by \\G\\ = Tr(GG'^f'^. Thus 
it is sufficient to consider only the pure states with 
standard Schmidt form given by Eq. It is easy 

to see that p = IV") ("01 = Y,i,j ^'^'^ 
p^^ = J2ij y/JWj\o-jbi) {a*bj \ ■ Then we arrive at 

j i 

= ||G®Gt|| = ||G|p = (^VMlf. (4) 

i 

where G = y^^ja*) {hj\. In this derivation we have 
used the unitarily invariant property of the trace norm 
when applying the elementary column transformation: 
{a*hj\ — > {hja*\ in the derivation of the first formula. 
The last formula is obtained from the observation that 

GGt - Y.^,,.J^W]\<){b^\ ■ \bj){a*\ = E.M.k*)(«:i, 

the property of the trace norm \\P ® Q\\ = \\P\\ ■ \\Q\\ 
and the fact that |1G|| is the sum of the square root of 
eigenvalues pi of GG'^ . 

Another complementary operational criterion for sep- 
arability called the realignment criterion is very strong 
in detecting many of BES 0, and even genuinely tri- 
partite entanglement 1^. Recently there has been con- 
siderable progress in the further analysis, and in finding 
stronger variants and multipartite generalizations for this 
criterion 0. We recall that this criterion states that a 
realigned version 7?,(p) of p should satisfy ||7?.(p)|| < 1 for 



any separable state p. TZ{p) is simply 7l{p)ij^ki = Pikji 
where i and j are the row and column indices for the sub- 
system A respectively, while k and I are such indices for 
the subsystem -B ll^(p)|| is also shown to be LU 

invariant in 0. One has Ti{p) = X^i j -y/MiAij joia*) {b*bj\ 
for the state Eq. (|2J), as follows easily from the definition. 
Similar to l@J one has 

mp)\\ = llE^/MIl«.)(&■l®Ev^l«*)(^J■Nl 

i j 

= ||G®G*|| = l|G!p = (EV7^)'- (5) 

i 

where G — J2iVl^\'^i) (b*\- The last formula fol- 
lows from the observation GG^ = j y/PiPj\ai) {b*\ ■ 

\bj) {a] \ = {ai\- 

We now derive the main result of this Letter. 

Theorem: For any m®n (to < n) mixed quantum state 
p, the concurrence C{p) satisfies 

C{p) > J . f (max(||p^^||,H7^(p)H) - l) . (6) 

Proof. — To obtain the desired lower bound, let us 
assume that one has already found an optimal de- 
composition 'Y^j^PiP^ for p to achieve the infimum of 
G{p), where are pure state density matrices. Then 
G{p) = Yl.PiGip') by definition. Noticing that Wp'^^W < 
E^P^\\ipT'\\ and ||7^(p)|| < E^P^mP')\\ duc to the 
convex property of the trace norm, one needs to show 

G{p') > ^2/{m{m-l)){\\{p')'^^\\ - 1) and G{p') > 
^2/ (m{m — l))(||7?.(p*)|| — 1). For a pure state one 

has \\n{p')\\ = WY^W = iJ2kVJ^? from Eqs. Q and 
(O , where y^/Ife are the Schmidt coefficients for the pure 
state p\ From the expression of Eq. Q it remains to 
prove that 

i<j ^ ' k 

= t ^ i, (EV7^f. (7) 

TO m — 1) ^-^ 

where we have used pi — 1. 

The verification of the inequality Eq. (J7|) is straightfor- 
ward: by summing over all of arithmetic mean inequali- 
ties piPj + pkPi > 2y/piPjPkPi for i < j and k < I, one 
gets 

'^^[PiP] + PktJ'i) > 2EE VM»A*jMfcMi 

i<j k<l i<j k<l 

= 2(E%/W^)'- (8) 

It is seen that the number of appearance times is m(TO— 1) 
for the term piPj on the Ihs of Eq. Therefore Eq. {Tj) 
is confirmed and the conclusion Eq. (|SJ| is proved. 
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The most prominent feature of the Theorem is that it 
allows to obtain an analytical lower bound for the con- 
currence without any numerical optimization procedure. 
The bound has the same range as C{p) and goes from 
to A/2(m — l)/m for pure product states and maximally 
entangled pure states, respectively. One can of course 
renormalizes the maximum value C{p) to be 1 with a 
change of the corresponding constant factor. 

We highlight some of the benefits of this new bound. 
Firstly, it serves to detect and gives a lower bound of 
concurrence for all entangled states of two qubits and 
qubit-qutrit system. This is so because the PPT cri- 
terion is necessary and sufficient for separability in the 
two cases 0. Secondly, it generalizes to bipartite sys- 
tems of arbitrary dimension a relation given in |20j 
which is only valid for two qubits case, that the con- 
currence is lower bounded by the negativity 0| (de- 
fined to be IIp"^-*!! — 1). Thirdly, for any qubit-qudit sys- 
tem our bound can contribute an analytical lower bound 
for EOF which is a convex function of the concurrence, 
E{\xP)) = i?2((l + ^1-C^m))l2) where 7?2(.) is the 
binary entropy function 11]. In fact our bound can fur- 
nish a lower bound of EOF E{p) for arbitrary bipartite 
state p Given any monotonously increasing, convex 
function £ satisfying £{C{\ip))) < —Yirp,r log2 /ir, one has 
E{p) > £{C{p)), with the rhs bounded from below by our 
bound Eq. ©. Next we consider some examples to illus- 
trate further the tightness and significance of our bound. 

Example 1: Isotropic states 

Isotropic states |2l|,|23| are a class oiU ®U* invariant 
mixed states in d x d systems 

P'^ - ^ (/- |vI/+)(*+|) +F|vl/+)(vl/+|, (9) 

where 1^-+) = VV^^Eti I") and F = 
satisfying Q < F < 1, is the fidelity oi pp and |\I'+). 
These states were shown to be separable for F < 1/d [21| . 
It is shown in [i[li that = ||7^(pF)|| = dF for 

F > 1/d. The concurrence C{p) for this class of states is 
recently derived in ^ to be y^2d/(d - 1){F- 1 /d) by an 
extremization procedure. An application of our Theorem 

gives C{p) > ^2/{d{d-l)){dF - 1) = C{p). Thus the 
bound gives surprisingly exact value of the concurrence 
for this sort of states. 

Remark: One can see that the equality of Eq. (O holds 
when 1-0) are product states or maximally entangled 
states (MES) (all are equal). Thus our bound will 
be tight if an optimal decomposition for achieving con- 
currence only involves product states and MES, and also 
attains the value of ||yO^^|| or ||7?.(p)||. Roughly speaking, 
the difference between our lower bound and the exact 
value of concurrence will be small if there are few de- 
viations from these two types of states in the optimal 
ensemble decomposition. Exact estimation of this differ- 
ence would be an interesting subject for future study. In 



the case of isotropic states, it is shown in |23| that the 
optimal decomposition falls exactly into this class and 
the concurrence is just our bound. 

Example ^."3x3 BES constructed from unextendible 
product bases (UPB) 

In l^^l, Bennett et al. introduced a 3 x 3 BES from 
the following bases: 
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from which the density matrix could be expressed as 



1 



(10) 



1=0 



A simple calculation gives \\p'^'^\\ = 1 and ||7?.(p)|| = 1.087 
0, therefore C{p) > 0.05 according to the Theorem. 
This shows that the state is entangled. 

When BES are constructed from the UPB given 

by = \vj] 



\V2j mod 5/ 



V, = 7V(cos(27rj75),sin(2^j75),/i), 



U = 0,...,4) with 
with j = 0,...,4, 

h = a/i + V5/2 and N ^ 2/^5 + VS, then the PPT 
state of Eq. ^ gives Wp'^-^W = 1 and ||7^(p)|| 1.098 
7], therefore C{p) > 0.056 according to the Theorem, 
which identifies this BES. 

It is conjectured by Audenaert et al. that the optimiza- 
tion method for concurrence is a necessary and sufficient 
for separability when one considers all possible complex 
linear combination of the concurrence- vectors 5j. Our 
numerical verification suggests a disproval for this con- 
jecture, because of the failure to identify entanglement 
by applying their optimization method for the above two 
UPB states. Thus the direct estimation method for con- 
currence m may not be able to detect all entangled 
states through numerical optimizations. Here our The- 
orem complements other existing approaches to make a 
quite good estimate of entanglement for BES. 
Example 3: Horodecki's 3x3 entangled state 
A mixed two qutrits is introduced in |25j : 



where 



a+ = -(|0)|1)(0|(1| + |1)|2)(1|(2|- 
a = i(|l)|0)(l|(0| + |2)|l)(2|(l|- 

|2)|2». 



|*+) = 



V3 



(|0>|0> 



11 



(11) 



|2)|0)(2|(0|), 
|0)|2)(0|(2|), 

(12) 



In 25] Horodecki et al. demonstrate that the states Eq. 
(|ll|l admit a simple characterization with respect to the 
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parameter 2 < a < 5: separable for 2 < a < 3; bound 
entangled for 3 < a < 4; free entangled for 4 < a < 5. 
It is computed by using the realignment criterion in Q 
that |i7^(cr„)|| = (19 + 2V3a2 - 15a+ 19)/21 and one 
can recognize all the entangled states for 3 < a < 5. 
One can obtain further that ||cr^'*|| = 1 for 2 < a < 4 
and |lcr;^-*|| = (2 + V4a2 _ 20a + 41)/7 for A < a < 
5. The refore one has Cjag) > l/V3(||7^(cra)|| - 1) = 
2-\/3(-\/3a^ — 15a + 19 — l)/63 due to the observation 
that ||7^(crQ)|| is always greater than ||crj^|| in the en- 
tangled region 3 < a < 5. 

However, the concurrence optimization procedure pro- 
posed in |l3| can only identify the entangled states for 
3.52 < a < 5 [l3- This suggests that Mintert et al.'s 
methods may not be necessary and sufficient for de- 
tecting entanglement. A rough comparison with the re- 
sult of (13j shows that our lower bound is much better 
than their optimized bound in the entangled region of 
3 < a ^ 4.75, though a little bit weaker than theirs in 
the region 4.75 ^ a < 5. 

We remark that, like any other known approaches, 
there are also some drawbacks for our estimation. Our 
lower bound cannot detect all the entangled states due 
to limitation of the PPT criterion and the realignment 
criterion. For example, it can neither recognize the 2x4 
Horodecki BES |2^ , which instead can be detected by the 
methods of 0|, nor give the exact value of concurrence 
for 2 qubits known from |3. Il3|. 

In summary, we have provided an entirely analytical 
formula for lower bound of concurrence, by making a 
novel connection with the known strong separability cri- 
teria. The bound leads to actual values of concurrence 
for some special class of quantum states. One only needs 
to calculate the trace norm of certain matrices, which 
avoids complicated optimization procedure over a large 
number of free parameters in numerical approaches. The 
formula also permits to furnish lower bounds of EOF for 
arbitrary bipartite quantum state. This complements the 
nice result of Wootters for 2 qubits, as well as a number 
of existing optimization methods for concurrence. Profit- 
ing from the strong realignment criterion, our bound can 
give easy entanglement evaluation for many BES, which 
fail to be recognized by the regular optimization meth- 
ods. This shows that our method can serve as a powerful 
tool for investigating both static and dynamical entangle- 
ment properties in realistic quantum computing devices. 
As applications the method could be used in indicating a 
possible quantum phase transition for condensed matter 
system, and in analyzing finite size or scaling behavior of 
entanglement in various interacting quantum many-body 
systems. 
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